Thin superconducting disk with 5-dependent Jc : Flux and current distributions 
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The critical state in a superconducting thin circular disk with an arbitrary magnetic field depen- 
dence of the critical sheet current, Jc{B), is analyzed. With an applied field Ba perpendicular to 
the disk, a set of coupled integral equations for the flux and current distributions is derived. The 
equations are solved numerically, and flux and current proflles are presented graphically for several 
commonly used Jc{B) dependences. It is shown that for small Ba the flux penetration depth can 
be described by an effective Bean model with a renormalized Jc entering the leading term. We 
argue that these results are qualitatively correct for thin superconductors of any shape. The results 
contrast the parallel geometry behavior, where at small Ba the -B-dependence of the critical current 
can be ignored. 



PACS numbers: 74.25.Ha, 74.76.Bz, 74.80. Bj05 
I. INTRODUCTION 

The critical state model (GSM) is widely accepted to 
be a powerful tool in the analysis of magnetic proper- 
ties of type-II superconductors. For decades there have 
been numerous theoretical works devoted to CSM cal- 
culations in the parallel geometry, i. e., a long sample 
placed in parallel applied magnetic field, Ba- More re- 
cently, much attention has also been paid to the CSM 
analysis of thin samples in perpendicular magnetic fields. 
For this so-called perpendicular geometry explicit ana- 
lytical expressions for flux and curreat, distributions have 
been, obtained for a long thin stripH'H and thin circular 
diskBij assuming a constant critical current (the Bean 
model). 

From experiments, however, it is well known that the 
critical current density jc usually depends strongly on 
the local flux density B. This dependence often hinders 
a precise interpretation of various measured quantities 
such as magnetization, complex ac susceptibility,QS and 
surface impedance.EI It is therefore essential to extend the 
CSM analysis to account for a i?-dependence of jc. 

In the parallel geometry extensive work has already 
been carried out, and exact-, rjesults for the flux det 
profiles and magnetization,E3"E3, as well as ac losses 
have been obtained for different jc (i?)-dependences. In 
the perpendicular geometry the magnetic behavior is 
known to be qualitatively different. In particular, due 
to a strong demagnetization, the field tends to diverge at 
the sample edges, and the flux penetration depth and 
ac tojses follow different power laws in Ba for small 
i3a.Ela Unfortunately, the theoretical treatment of the 
perpendicular geometry is very complicated, and we are 
not aware of any explicit expressions obtained for the 
CSM with a i?-dependent jc- However, it is possible to 
derive integral equations relating the flux and current 
distributions .□ Such equations have so far been obtained 
and s olvei d numerically only for the case of a long thin 

strip.ye 



In this paper, we derive a CSM solution for a thin 
circular disk characterized by an arbitrary jc{B). The 
solution is presented as a set of integral equations which 
we solve numerically. In this way we obtain the field 
and current density distributions in various magnetized 
states. We present results for several commonly used 
functions jc{B). A special attention is paid to the low- 
field asymptotic behavior. 

The paper is organized as follows. In Sec. Il^ the basic 
equations for the disk problem are derived. We consider 
here all states during a complete cycle of the applied field, 
including the virgin branch. Sec. Ill contains our numer- 



ical results for fiux and current distributions as well as 
for the fiux front position. A discussion of the results is 
presented. Finally, Sec. IV presents the conclusions. 



II. BASIC EQUATIONS 

Consider a thin superconducting disk of radius R and 
thickness d, where d <^ R, see Fig. 1^. We assume either 
that d > A, where A is the London penetration depth, 
or, if d < A, that }? /d <C R- In the latter case the 
quantity \\/-d plays a role of two-dimensional penetra- 
tion depth.E3 We put the origin of the reference frame 
at the disk center and direct the z-axis perpendicularly 
to the disk plane. The external magnetic field Bq is ap- 
plied along the z-axis, the z-component of the field in the 
plane z ~ Q being denoted as B. The current flows in 
the aziniuthal direction, with a sheet current denoted as 
•^(^) — /-d/2-^('"' '^) ^^' ■^here j is the current density. 

To obtain expressions for the current and flux dis- 
tribution we follow a procedure originally suggested in 
Ref. g and then generalized in Ref. 14 for the case of in- 
dependent Jc in a thin strip. The procedure makes use 
of the Meissner state distributions for B and J. 

In the Meissner state, where B,,^ inside the disk, the 
field outside the disk is given byH'Q 



BM{r,R)^Ba 



2Ba 



R . [R 

= — arcsm — 

Vr2 - i?2 \r 



and the current is distributed according to 

4Ba r 



JM(r,i?) 



TTAiQ ^ija _ J.2 



r < R 



, (1) 



(2) 



Be — ^0 Jco/2 , J^cO = ^c(5 — 0). 



(8) 



Note that due to the similar form of the function Jm (?", 
in (0) and the Meissner-state current in the strip cas^ 
our weight function G{r, Bg) is also similar to that for a 
strip, see Ref. n4 

From Eqs. (0) and (M) it then follows that the current 
distribution in a disk is given by 




\ 



FIG. 1. Schematic of a thin circular disk in a perpendicular 
applied field. The current flows in the azimuthal direction. 



A. Increasing field 

We begin with a situation where the external field Ba 
is applied to a zero-field-cooled disk. The disk then con- 
sists of an inner flux-free region, r < a , and of an outer 
region, a < r < R, penetrated by magnetic flux. Accord- 
ing to the critical state model, the penetrated part will 
carry the critical sheet current Jc corresponding to the 
local value of magnetic field. 



Jir) 



-Jc[B{r)], a<r < R 



(3) 



Now following Refs. ^,|lj, we express the field and current 
as superpositions of the Meissner-state distributions, ffl) 
and (|), i.e.: 



^min(r,i?) 



B{r)= / dr'BM{r,r')G{r',Ba). 

J a 



Jir) 



niax(a,r) 



dr'JM{r,r')Gir',Ba), 



(4) 



(5) 



where G{r,Ba) is a weight function. Since B{r) and 
Bm (r, i?) —> i?a at r ^ oo we have the normalization 
condition 



f drG{r,Ba) = l. 

J a 



(6) 



Substituting Eq. (|) and Eq. (|) into Eq. (|) yields an 
integral equation for the function G{r, Ba) which can be 
inverted to obtain 



-^-^'^-kil' 



dr' Jc[B{r')] 



Vf 



(7) 



J{r) 



2r 

-Jc\B(r) 



dr' 



'MB{r')] 



r < a 
a <r < R 



(9) 



This equation is supplemented by the Biot-Savart law, 
which for a disk reads, a 



B(r) = B„ 



27r 



F{r,r')J{r')dr 



(10) 



Here F{r,r') = K{k)/{r + r') - E{k)/{r - 

r'), where k{r,r') = 2-\/rr'/(r + r'), while K 

and E are complete elliptic integrals defined as 

E{z) = J^^'^[l~z^cos^x)]^^^dx and K{z) = 



-1/2 



dx. 



j;^'[l-z^cosHx)] 

The relation between the flux front location a and ap- 
plied field Ba is obtained by substituting Eq. (0) into 
Eq. (i), giving 



Ba = Br 



dr' Jc[B{r')] 



a V^ 



J. 



cO 



(11) 



For a given Ba and for a specified Jc{B) we need to 
solve the set of three coupled equations (p|)-(pd]). In the 
case of B-independent J^ the Eq. (|l^) acqrures the sim- 
ple form 



a/R^l/cosh{Ba/B,), 



(12) 



and the Eqs. (S) and (H) lead to the Bean-model results 
derived in Refs. ^ and^ 

Note that the equations can be significantly simplified 
at large external field where a — > proportionally to 
ex.p{—Ba/Bc)- Then B{r) is determined by the single 
equation 



Bir)^Ba~^J Fir,r')MB{r')]dr' ., 



(13) 



where 



following from Eq. (hd 



B. Subsequent field descent 

Consider now the behavior of the disk as Ba is reduced 
after being first raised to some maximal value Bam- Let 
us denote the flux front position, the current density and 
the field distribution at the maximum field as Om, Jmir) 
and Bmir), respectively. Obviously, Jmif), Bm{r), and 
am satisfy Eqs. (i)-P). 




FIG. 2. Flux density profile during the applied field descent 
from a maximum value Bam- 

During the field descent from Bam the flux density be- 
comes reduced in the outer annular region a < r < R, 
see Fig. g. The central part of the disk r < a remains 
frozen in the state with Ba = Bam- Let us specify the 
field and current distributions in this remagnetized state 
as 

Bir) = Bmir) + B{r), J{r) = Jm{r) + J{r), (14) 

and derive the relation between B{r) and J{r). For that 
one can use a procedure similar to the one described 
in Sec. [I A. The only difference is that in the region 
a < r < R we now have to use J{r) = +Jc[B{r)]. In 
this way we obtain 



Mr), 



a <r < R, 



J{r) 
where we define 

Jc{r) = Jc[Bm{r) + B[r)] + Jc[Bm{r 



(15) 



Note that the function Jdr) depends on the coordinate 
only through the field distributions Bm{r) and B{r). In- 
stead of Eq. (0), the additional current satisfies 



Jir) 



with the complementary equation 



T^ Ja 


Vr'2 


-r2 J,{r') 
- a2 r'2 - r2 ' 


r < a 


Mr), 






a <r < R 

(16) 



B{r) = Ba 



Ba 



Mo 

27r 



F{r,r')J{r')dr' . (17) 



Furthermore, similarly to Eq. (pT|), we have 



Ba — Ba 

bT 



dr 



\fr^ 



Mr) 

JcO 



(18) 



which completes the set of equations describing the re- 
magnetized state. Again, for ^-independent Jc the equa- 
tions reproduce the Bean-model resultaJ'B 

If the field is decreased below —Bam the memory of 
the state at Ba = Bam is completely erased, and the so- 
lution becomes equivalent to the virgin penetration case. 
If the difference Bam — Ba is sufficiently large then a ^ 
rapidly, and the critical state J{r) = Jdr) is established 
throughout the disk. In this case the field descent is de- 
scribed by Eq. (O) with the opposite sign in front of the 
integral. 

We emphasize that the expressions derived here (Sec. II 
A and B) are readily converted to the long thin strip case. 
This is due to the similarity of (0) and the expression for 
the Meissner-current in a strip. 



J{x) 



2Ba 



Mo \puP^—x^ 



-w < X < w, 



(19) 



where x is the coordinate across the strip. Thus, making 
in this paper the substitutions 



X, R ^ w, F{r, r') 



2x' 



J2^ ^'= 



fJ-oJ, 



cO 



one immediately arrives at the set of equations valid for 
a thin long strip. In that case some of the integrals can 
be done analytically to yield the expressions obtained in 
Ref . |l|. 

A difference between the derivation in Ref. 113 and the 
present one is that for decreasing fields we calculate only 
the additional field B{r) rather than the total field B{r). 
This allows us to use only one weight function (0) to cal- 
culate the flux distributions both for increasing and de- 
creasing fields. This simplifies the numerical calculations 
significantly. 



C. Numerical procedure 

Given the Jc(i?)-dependence, the magnetic behav- 
ior is found by solving the derived integral equations 
(m)"(|11[) numerically using the following iteration pro- 
cedure. With Ba increasing a flux front position a is 
first specified, and an initial approximation for B{r), 
e.g., the Bean model solution, is chosen. At each step 
the nth approximation, i3'^")(r), is used to calculate 
j(»)(r) from Eq. (|) and si"^ from Eq. (0). They are 
then substituted into Eq. (Iln) yielding the next approx- 
imation, i3'^"+^^(r). The iterations are stopped when 



(R-^Jdr [B("+i)(r)-B(")(r)]^y ^ < IQ-^B^. With 

Ba decreasing, the same procedure is used to find first 
Jmir), Bmir), Bam. for a given am- Then, Eqs. @-(|l|) 
are solved for a fixed a yielding the functions B{r) and 
J{r) and also the applied field Ba- 



III. FLUX AND CURRENT DISTRIBUTION 

A. General features 

In the numerical calculations we used the following de- 
pendences Jc{B), 
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FIG. 3. Flux density (a) and current (b) profiles for flux 
penetration into a virgin state. Solid lines are calculated for 
the Bean model, while other ones represent the Kim model, 
Eq. (|2C|), with different Jco and Bo- These parameters are 
chosen so that the flux front position a = 0.2i? and applied 
fleld Ba are the same for all the proflles. The corresponding 
Jc (S)-dependences are shown in the panel (c). The current 
is normalized to the Bean model critical current J?''^". 



Jc = Jco/(l + |B|/So) 

Jc = Jeo exp(-|B|/Bo) 



(Kim model), (20) 

(exponential model). (21) 



Shown in Fig. (a,b) are the field and current distri- 
butions for increasing field with ajB = 0.2 for the Kim 
model with different parameters Jco and Bq. They are 
chosen in a way to keep the position a fixed for all the 
curves for a given value of Ba- This allows us to follow 
the variations in the profile shape as the S-dependence 
of Jc changes. The chosen parameters Jco and Bq cor- 
respond to the set of Jc(i3)-curves shown in Fig. ^ (c). 
The Bean model results are also plotted in Fig. 0. 
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FIG. 4. Graphs similar to Fig. H, but here for the expo- 
nential model, Eq. (pi[). 

Several major deviations between the Kim and the 
Bean model can be noticed. In the Kim model we see 
that (i) the current J(r) is not uniform at a < r < i? - 
it is minimal at the disk edge where |i?| is maximal; (ii) 
the current has a cusp-like maximum at r = a since the 
magnetic field vanishes at this point with infinite deriva- 
tive; (in) compared to the Bean model, the B(r) profiles 



are steeper near the flux front, whereas the peaks at the 
edges are less sharp. 
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FIG. 5. Flux density and current profiles for decreasing 
Bo,- The graphs show the Kim model with Bo/ Be ~ 3 
and the Bean model for Ba = Bam ~ 2ABc (solid), 
Ba = 1.3Bc (dashed), Ba = 1.3Bc (dotted), and Ba = -2ABc 
(dash-dotted). 



Qualitatively similar results are obtained for the ex- 
ponential model, see Fig. ^. Also here, by changing the 
model parameters one can produce a variety of flux and 
current profiles which are quite different from the Bean 
model predictions. When comparing the Kim and expo- 
nential model, however, it turns out to be very difficult 
to find clear distinctions. 

During field descent the B- and J-profiles become more 
complicated. For brevity we show only profiles for the 
Kim model with Bq/Bc = 3 and for the Bean model at 
different values of Ba, see Fig. ||. Again, the Kim model 
gives a non-uniform current density at r > a. Contrary 
to the increasing-field states, the current density can now 
either decrease, or increase towards the edge depending 
on Ba- 

Figure ra shows profiles for fully-penetrated decreasing- 
field state. In the Bean model the current remains con- 
stant, while the profile of the flux distribution is fixed, 
although shifted according to the applied field. In con- 
trast, the Kim- model profiles are strongly dependent on 
Ba- There is a peak in the current profile and an en- 
hanced gradient of B(r) near the point where B = 0. 




r/R 

FIG. 6. Flux and current density profiles during field de- 
scent after first increasing Ba to a very large Bam- It is as- 
sumed that critical state is established throughout the disk. 
Solid lines show results for the Kim model with Bq — Be 
and Ba/Bc = -1.4 (curve 1), Ba/B^ = -0.7 (curve 2), and 
Ba ~ (curve 3). The dotted lines show the Bean model 
profiles, their shape being independent of Ba- In contrast, 
the shape of the Kim model profiles depends strongly on Ba- 



B. Flux penetration depth 



To analyze quantitatively the role of a _B-dependent 
Jc let us consider the position of the flux front during 
increasing field. A circle with radius a then limits the 
Meissner region i? = 0, and is also the location of max- 
imum gradient in B. These features can be measured 
directly in experiments on visualization of nmgnetic flux 
distribution, e. g., magneto-optical imaging.cj 

Let us first recall the CSM expression for the flux fjmnt 
location in a long circular cylinder in a parallel ficldEj, 



a 
R 



1 



^o-R 



■ dB 



(22) 



At small applied fields, Ba, it can be expanded as 



a Ba fJ-oR .,, . f Ba 



(23) 



where for a cylinder Be = l^ojcoR- Note that the B- 
dependence of jc enters the expansion first in the second- 
order term. Consequently, for a long cylinder the low- 
field behavior of a is well described by the Bean model 
where the penetration depth increases linearly with the 
applied field. 

For a thin disk, the penetration of flux proceeds differ- 
ently. In the Bean model the location of the flux front, 
Eq. ([l^), is for small Ba given by 



(24) 



a/i? « 1 - (1/2) (Ba/Bc)' 



For an arbitrary i?-dependence of Jc the expression (|l] 
relating a and Ba cannot easily be expanded in powers of 
the ratio Ba/Bc- The physical reason for this is the sin- 
gular behavior of the magnetic field near the disk edge. 
There, the local field diverges at any finite Ba, and an 
expansion of Jc{B) in powers of B is not everywhere con- 
vergent. To clarify the behavior of the fiux front we have 
therefore performed numerical calculations of the depen- 
dences a{Ba,Bo). Shown in Fig. are the results for 
the Kim (upper panel) and the exponential (lower panel) 
models. Note that the limit of large Bq represents the 
Bean model. 

For small Ba all the models seem to yield a parabolic 
relation between the penetration depth and the applied 
field. This is illustrated in more detail in Fig. 0, where 
all the graphs in the log-log plot have a slope of 2 in the 
low-field region. We therefore conjecture that any B- 
dcpendcnce of Jc leads to the same quadratic law ( 24 ) as 
for the Bean model, although with different coefficients 
in front oi[Ba/Bcf. 




FIG. 7. Flux front location a as a function of increas- 
ing applied field, Ba, and the charcteristic field, Bq. The 
calculated results are presented for the Kim model (upper 
panel), Eq. (|20|), and for the exponential model (lower panel), 
Eq. @). 
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FIG. 8. Reduced fiux penetration depth, 1 — ajR, in the 
virgin flux-penetrated state versus the applied field, Ba- For 
all three models the dependences have quadratic asymptotic 
behavior at low fields. 

The overall behavior of the penetration depth can be 
fitted well by the full form Eq. ( [l2[) , provided one makes 
the substitution Be -^ B^^ , i-e., 

a/i? = l/cosh(Ba/Bf) . (25) 

Wc find that the effective Be satisfies the relation 



Bf 
Be 



1 



Be 
Bo' 



(26) 



if the ratio Bq/Bc is of the order 1, or larger. Here 
a — 0.42 for the exponential model, and a = 0.36 for the 



Kim model. The same relations ( |25D and (26) are found 
to hold true for a long thin strip with a = 0.60 and 0.51 
for the exponential and the Kim model, respectively. 

We believe that for many purposes a Bean-model de- 
scription with an effective critical current is appropriate 
for thin samples of any shape, both in applied field and 
under transport current. Indeed, strong demagnetization 
effects always lead to a divergence of magnetic field at the 
sample edge. This implies that in the sample there is al- 
ways present a wide range of B values up to infinity. As 
a result, the sample behavior is determined by the whole 
Jc (i3)-dependence. In particular, the value Jc(0) is not 
governing the magnetic behavior ot-thin samples, even 
when the applied field is very small.ll3 



IV. CONCLUSION 

A set of integral equations for the magnetic flux and 
current distributions in a thin disk placed in a perpen- 
dicular applied field is derived within the critical state 
model. The solution is valid for any field-dependent crit- 
ical current, Jc{B). By solving these equations numer- 
ically it is demonstrated that both the flux density and 
current profiles are sensitive to the Jc(B)-dependence. 
In particular, compared to the Bean model, the B{r)- 
profiles are steeper near the flux front, whereas the peaks 
at the edges are less sharp. 

Since the local magnetic field at the disk edge is di- 
vergent for any value of the applied field, Ba, a field 
dependence of Jc affects the flux distribution even in the 
limit of low Ba- Our numerical calculations show that 
the flux penetration depth at small fields has the same 
quadratic dependence on Ba as for the Bean model, how- 
ever with different coefficient. The overall behavior of 
the ffux penetration depth is well described by the Bean 
model expression with an effective value of the critical 
current. These results are believed to be qualitatively 
correct for thin superconductors of any shape. The be- 
havior differs strongly from the case of a long cylinder 
in a parallel held, where the front position at low Ba is 
not affected by the B-dependence of the critical current 
density. 



ACKNOWLEDGMENTS 



The financial support from the Research Council of 
Norway is gratefully acknowledged. 



^ E. H. Brandt, and M. Indenbom, Phys. Rev. B 48, 12893 

(1993). 
^ E. Zeldov, J. R. Clem, M. McElfresh, and M. Darwin, Phys. 

Rev. B 49, 9802 (1994). 
^ P. N. Mikheenko and Yu. E. Kuzovlev, Physica C 204, 229 

(1993). 
^ J. Zhu, J. Mester, J. Lockhart, and J. Turneaure, Physica 

C 212, 216 (1993). 
^ J. R. Clem and A. Sanchez, Phys. Rev. B 50, 9355 (1994). 
^ Some final expressions obtained in Ref. ra are not correct, 

see discussion in Ref. H. 
■^F. Gomory, Supercond. Sci. Technol. 10, 523 (1997). 
B. J. Jonsson, K. V. Rao, S. H. Yun, and U. O. Karlsson, 
Phys. Rev. B 58, 5862 (1998). 
® B. A. Willemsen, J. S. Derov, and S. Sridhar, Phys. Rev. 

B 56, 11989 (1997). 
^° J. R. Clem, J. Appl. Phys. 50, 3518 (1979). 
" M. Forsthuber and G. Hilscher, Phys. Rev. B 45, 7996 

(1992). 
^^ M. Xu, D. Shi, R. F. Fox, Phys. Rev. B 42, 10773 (1990). 
" T. H. Johansen and H. Bratsberg, J. Appl. Phys. 77, 3945 

(1995). 
^^ J. McDonald and J. R. Clem, Phys. Rev. B 53, 8643 (1996). 
^* An alternative approach tp-|1pEat the problem is to carry 
out flux creep simulations.tJO Here the GSM results can 
be found as a limiting case of a large creep exponent. 
^^ In samples of finite thickness d the field divergency at 
the sample edge is cut off at values of the order of 
-Be ln(w/d). Thus, the low-field behavior is determined by 
Jc(-B)-dependence at fields less than this value. 
" E. H. Brandt, Phys. Rev. B 55, 14513 (1997). 
^* E. H. Brandt, Phys. Rev. B 58, 6506 (1998). 
^^ P. G. de Gennes, Superconductivity of Metals and Alloys 

(Benjamin, New York, 1966). 
'^° Anjah B. Riise, T. H. Johansen, H. Bratsberg, M. R. Kob- 
lischka, and Y. Q. Shen, submitted to Phys. Rev. B.; M. R. 
Koblischka and R. J. Wijngaarden, Supercond. Sci. Tech- 
nol. 8, 199 (1995). 



